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Abstract 

Background: HIV-positive women are at increased risk of human papillomavirus (HPV) infection, and, especially, 
multiple infections compared to HIV-negative women. Whether certain HPV types have a tendency to cluster in 
multiple infections beyond or below what would be expected by shared risk factors (e.g., sexual behavior and the 
degree of immunosuppression) is unclear. We, therefore, investigated clustering patterns of 44 HPV types in 
HIV-positive women from Kenya. 

Findings: HPV status was assessed on cervical scrapings from 498 women using GP5+/6+ PGR and reverse line 
blot. Logistic regression was used to model type-specific HPV positivity, adjusted for age, specific HPV type prevalence, 
CD4, combination antiretroviral therapy, and, in the Full Model, individual-level random effects that represent 
unobservable risk factors common to all HPV types. We found a modest excess of women with co-infections with 2 
HPV types (1.12; 95% credible intervals: 1.03-1.21) in the Full Model but no significant associations of individual types. 
No significant deviations of observed/expected counts were observed for any 2-way combination of HPV types at either 
the chosen level of significance, p = 0.00005, or at p = 0.01 . Findings were substantially similar when women with CIN2/3 
were excluded and when they were stratified by use of anti-retroviral therapy or CD4 count. 

Conclusions: HPV co-infections occurred at random in the cervix of HIV-positive women as previously found in 
HIV-negative women. The removal of HPV types through vaccination should not result, therefore, in an increase or 
decrease in the prevalence of HPV types not targeted by vaccination in immunosuppressed women. 
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Findings 

Introduction 

HIV-positive women are at increased risk of infections 
with human papillomavirus (HPV) compared to HIV- 
negative women [1,2]. Although multiple HPV infections 
are particularly common in HIV-positive women, it is 
not clear whether certain HPV types have the tendency 
to cluster in multiple infections beyond or below what 
would be expected by shared risk factors, e.g., sexual be- 
havior and the degree of immunosuppression. This 
could be relevant for the evaluation of the effects of 
HPV prophylactic vaccines. The removal of certain HPV 
types through vaccination could, in theory, indirectly 
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increase or decrease the prevalence of other untargeted 
types. 

Previous studies among HIV-negative women [3-7] 
and men [7,8] have reported a general tendency of HPV 
types to cluster in multiple infections. However, there 
has been no evidence for specific HPV types to be found 
together more or less often than expected by chance 
with other types. Significant excesses of a few HPV type 
combinations have been observed, but they have been 
demonstrated to derive from diagnostic artefacts, e.g., 
cross -hybridization of closely homologous HPV types 
using enzyme immunoassay for genotyping [8] or tech- 
nical limitations of the HPV detection methods used, e.g., 
indirect measure of HPV52 by using a mixed probe using 
Roche Linear Array [3]. 

To date, the clustering patterns of multiple HPV infec- 
tions in HIV-positive women have been evaluated in only 
one study based in the United States [9]. Chaturvedi et al. 
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[9], however, only assessed clustering at the phylogenetic 
clade level, not individual HPV type level We, therefore, 
have investigated this issue in a cross-sectional study of 
HIV-positive women in Nairobi, Kenya. 

Materials and methods 

The study in Kenya has been described in previous publi- 
cations [10-12]. Briefly, in 2009, 500 HIV-positive women, 
aged 18-55 years, were enrolled with the primary aim 
of comparing cervical cancer screening methods. After 
obtaining a written informed consent, cervical exfoliated 
cells were successfully collected from 498 women to per- 
form HPV testing, and HPV DNA was assessed using a 
general GP5+/6+ primer-mediated PGR to detect 44 HPV 
types with reverse line blot analysis for genotyping [13]. A 
venous blood sample was taken to measure CD4 count. 
Type-specific prevalence was described in detail [10]. HPV 
prevalence is slightly different in the present report [10], 
because six women who were HPV-positive only according 
to HPV generic mix probes were excluded. The study 
protocol was approved by the Ethical Review Committees 
of the Kenyatta National Hospital, Kenya; the University 
of Washington, USA; and the International Agency for 
Research on Cancer, France. 

Statistical methods used in the present report have 
been described in previous publications [3,5-8,14,15]. 
Briefly, a multivariate logistic regression was used to 
model type-specific HPV positivity. None of the sexual 
behavior variables collected by the enrollment question- 
naire was included in the models because they were not 
found to be significantly associated with HPV-positivity 
in a previous analysis [10]. 

Three models were used. The Basic model included 
age and specific HPV type prevalence only as covariates. 
The Adjusted model added CD4 count (as a continuous 
variable), use of combination antiretroviral therapy (cART) 
(never, <2 years, >2 years) and an interaction term between 
CD4 counts and cART use. This model adjusted for the 
fact that low CD4 counts were found to be associated with 
higher HPV positivity and that this association was weaker 
among women with a prolonged cART use [10]. The Full 
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Figure 1 Observed versus Expected occurrence for 2-way 
combinations of 44 HPV types, HIV-positive women in Kenya. 

Plus signs represent occurrences of HPV pairs. HPV pairs located in 
the upper triangle indicate positive clustering, while those located in 
the lower triangle represent negative clustering between the HPV 
types involved. None of the p-values for joint HPV infections was 
significant at the chosen significance level of 0.00005. No significant 
associations emerged even when considering a less conservative 
level, p<0.01. 
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model added individual-level random effects. Random ef- 
fects quantify and account for the fact that women with 
several concurrent HPV types may have a higher level of 
risk for HPV infection compared to women with few or 
no HPV types, and that type-specific HPV measurements 
in the same woman are correlated with each other. 
Individual-level random effects represent unobservable risk 
factors common to all HPV types [3]. 

A Bayesian approach with Markov Chain Monte Carlo 
simulation was used; estimates were reported as pos- 
terior means and 95% credible intervals (95% CI). 
Discrepancies between the data (observed counts of 
co-infections per subject) and the model (expected 
counts of co-infections per subject) were assessed by 
posterior predictive two-sided p-values and measured 



Table 1 Observed (O) to Expected (E) ratio of multiple infections with 44 HPV types, according to various models 









Basic model 




Adjusted model 


Full model 




No of HPV types 


0 


% 




0/E (95% CI)'' 




0/E (95% Cl)^ 


E^ 


0/E (95% Cl)^ 


0 


162 


32.5 


122.8 


1.32 (1.20-1.46) 


131.5 


1.23 (1.11-1.35) 


163.1 


0.99 (0.90-1.09) 


1 


145 


29.1 


175.2 


0.83 (0.81-0.85) 


168.7 


0.86 (0.84-0.89) 


151.2 


0.96 (0.92-1.00) 


2 


103 


20.7 


121.4 


0.85 (0.81-0.89) 


114.0 


0.90 (0.86-0.96) 


92.3 


1.12 (1.03-1.21) 


3+ 


88 


17.7 


78.6 


1.13 (0.98-1.29) 


82.8 


1.07 (0.93-1.22) 


90.4 


0.98 (0.87-1.11) 



O = Observed, E = Expected, CI = Credible Interval. 
^Controlling for age, specific HPV prevalence. 

^Controlling for age, specific HPV prevalence, CD4 counts, cART and interaction term between CD4 and cART. 
^As^ plus individual random effects. 
3+ means 3 HPV types or more. 



Vaccarella et al. Infectious Agents and Cancer 2013, 8:50 
http://www.infectagentscancer.conn/content/8/1/50 



Page 3 of 4 



by an observed-to-expected (O/E) ratio for each HPV 
co-infection. Since all possible two-way interactions 
between the 44 HPV types were included, this generated 
946 (44 X 43/2) statistical comparisons. To minimize errors 
due to multiple comparisons, the Bonferroni correction 
was used to set the p-values thresholds in order to assess 
statistical significance. With 946 multiple comparisons, the 
Bonferroni corrected p-value threshold for type-type asso- 
ciations was 0.05/946 = 0.00005. 

Results 

Overall HPV positivity in the 498 HIV-positive women 
was 67.5% (N = 336). Multiple infections were found in 
56.8% (N=191) of HPV-positive women (Table 1). 
According to the Basic model, there was an excess of 
HPV-negative women [O/E, ratio = 1.32; 95% Credible 
Interval (CI): 1.20-1.46]. O/E ratios for infection with 2 
and >3 HPV types were 0.85 (95% CI: 0.81-0.89) and 
1.13 (95% CI: 0.98-1.29), respectively. With the Adjusted 
model, the O/E ratio was 1.23 (95% CI: 1.11-1.35) for 0 
HPV types, 0.90 (95% CI: 0.86-0.96) for 2 HPV types, 
and 1.07 (95% CI: 0.93-1.22) for >3 HPV types. Inclusion 
of CD4 counts, cART use and their interaction in the 
model slightly reduced the discrepancy between ob- 
served and expected counts. When the Full model was 
used, no substantial difference was found between ob- 
served and expected counts, with the exception of a 
small but significant excess of double HPV infections 
(O/E ratio = 1.12, 95% CI: 1.03-1.21). 

No significant deviations of O/E counts were observed 
for each 2- way combination of 44 HPV types (Figure 1). 
None of the pairs of HPV types reached the chosen level 
of significance, p = 0.00005. No significant excess or def- 
icit of 2-way combinations of types was found when we 
used a less conservative level of significance (p = 0.01) 
(data not shown). The results of the analyses were sub- 
stantially similar when women with CIN2/3 were ex- 
cluded from the analyses and when stratified by cART 
use (never; ever) and CD4 count (< or >350 cells/(iL) 
(data not shown). 

Discussion 

The present study of HIV-positive women in Kenya found 
that there was no evidence of either an excess or a lack of 
clusters between specific HPV types, even though a slight 
excess of infections with 2 HPV types was observed even 
after controlling for all sources of common correlation be- 
tween HPV types (random effects). 

Our study is the first to evaluate clustering patterns 
of specific HPV types among HIV-positive women in 
Africa. The assessment of the associations between all 
the different HPV genotypes requires studies with a high 
statistical power. The present study was relatively small 
compared to previous studies in the general population 



[3,5,6,8]. Small size was, however, partly compensated by 
the larger number of multiple HPV infections in HIV- 
positive women. 

A strength of our analysis was the inclusion into the lo- 
gistic regression model of random effects at the individual 
level that allowed us to account for the unobservable risk 
factors shared by all HPV types. Non type-specific cluster- 
ing of HPV types has often been observed in previous 
studies of HIV-negative women prior to the inclusion of 
random effects. Non type-specific clustering has generally 
been explained by a correlation with high-risk sexual be- 
havior as all HPV types are sexually transmitted [3,5,6,8]. 
In our study of HIV-positive women from Kenya, however, 
multiple HPV infections were associated with the degree 
of immunosuppression and not with sexual behavior [10]. 
Immunosuppression diminishes the probability of HPV 
clearance and increases the risk of HPV reactivation [16]. 
Random effects, however, contributed substantially to 
reducing the spurious discrepancy between observed and 
expected numbers of infections, with the possible excep- 
tion of double infections. 

In conclusion, our present findings confirm in HIV- 
positive women previous results observed in HIV-negative 
women, providing further evidence on the lack of cluster- 
ing between carcinogenic HPV types in either immuno- 
competent or immunosuppressed women. 
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